Theoretical insight of reactive oxygen species scavenging mechanism in lignin waste depolymerization products

Apart from natural products and synthesis, phenolic compounds can be produced from the depolymerization of lignin, a major waste in biofuel and paper production. This process yields a plethora of aryl propanoid phenolic derivatives with broad biological activities, especially antioxidant properties. Due to its versatility, our study focuses on investigating the antioxidant mechanisms of several phenolic compounds obtained from renewable and abundant resources, namely, syringol (Hs), 4-allylsyringol (HAs), 4-propenylsyringol (HPns), and 4-propylsyringol (HPs). Employing the density functional theory (DFT) approach in conjunction with the QM-ORSA protocol, we aim to explore the reactivity of these compounds in neutralizing hydroperoxyl radicals in physiological and non-polar media. Kinetic and thermodynamic parameter calculations on the antioxidant activity of these compounds were also included in this study. Additionally, our research utilizes the activation strain model (ASM) for the first time to explain the reactivity of the HT and RAF mechanisms in the peroxyl radical scavenging process. It is predicted that HPs has the best rate constant in both media (1.13 × 108 M−1 s−1 and 1.75 × 108 M−1 s−1, respectively). Through ASM analysis, it is observed that the increase in the interaction energy due to the formation of intermolecular hydrogen bonds during the reaction is an important feature for accelerating the hydrogen transfer process. Furthermore, by examining the physicochemical and toxicity parameters, only Hs is not suitable for further investigation as a therapeutic agent because of potential toxicity and mutagenicity. However, overall, all compounds are considered potent HOO˙ scavengers in lipid-rich environments compared to previously studied antioxidants.


Introduction
7][8] Consequently, there has been a focus on developing various lignin-based materials, including bioplastics, supercapacitors, and batteries, to harness the technical lignin derived from pre-treatment steps. 9,10Lignin depolymerization is a crucial step for maximizing the economic value of the valorization process of lignin into value-added products. 11The depolymerization process involves the cleavage of b-aryl ether linkages, followed by subsequent reactions that give rise to several aromatic compounds, which are phenol, guaiacol, and syringol derivatives, that exhibit broadspectrum application in human health and industry. 12,135][16][17] In addition, many of them have been proven to possess many biological activities such as antifungal, antibacterial, and antioxidant. 18,19The versatile antioxidant nature of these compounds enables their potential application in the pharmaceutical industry and other sectors.It could be employed as a preservative in dried cereal to prevent polyunsaturated fatty acid oxidation. 20,21][24][25] From three major structures of lignin units, syringol is reported to have potent antioxidant activity comparable with vitamin C and E. 26 Furthermore, its derivatives exhibit better radical scavenging capabilities compared to other groups, as demonstrated through in vitro or in silico analyses. 27,280][31][32] Among other derivatives, 4allysyringol and 4-propenylsyringol are two of the many derivatives with under-explored antioxidant activities.Furthermore, these are included as staple derivatives of syringol produced from lignin depolymerization.4][35][36] Besides that, the presence of an unsaturated side chain is the main reason these compounds are selected for this work.A previous study claimed that it could increase the antioxidant potential of phenolic derivatives more than the parent structure. 37Moreover, the conjugative effect on this group results in good peroxyl radical scavenging activity in the lipid environment. 29,38This could lead to the development of neuroprotective drugs in the future. 39owever, a deep investigation of the antioxidant mechanism of these two compounds is still not conducted.Thus, this work mainly consists of elucidating the radical scavenging of these compounds using computational studies.Along with 4-allylsyringol (HAs) and 4-propenylsyringol (HPns), syringol (Hs) and 4-propylsyringol (HPs) are also investigated to serve as the comparison structures (Fig. 1).We choose the hydroperoxyl radical (HOOc) in this study as a target of scavenging for several considerations.This radical is the simplest form of the peroxyl radical (ROOc) group, which is involved in oxidative damage, such as lipid peroxidation and brain tumor build-up. 40,41In addition, HOOc is basically the protonated form of the superoxide radical and is proposed to be more reactive than its conjugated base in living systems despite its low availability. 42owever, the main reason is its mild reactivity, which can be seen from its half-life value. 43Thus, it can give a clear difference in the antioxidant efficiency of the tested compounds.Meanwhile, highly reactive radicals, such as hydroxyl radicals, tend to give diffusion-limited rates for any compounds, which biases the result.Thus, by modeling the scavenging mechanism of syringol derivatives with HOOc, we will provide reliable kinetics data to gain insight into the structure-activity relationship for designing potent antioxidant compounds in specic environments.
In scavenging free radicals, antioxidant compounds provide various reaction mechanisms with different reactivity.From our previous research, one observed phenomenon is the predominance of the HT (hydrogen transfer) mechanism by phenolic groups in antioxidant compound. 44Meanwhile, other fragments, such as benzyl groups, exhibit higher activation energy.
Thermodynamic and kinetic assessments are still insufficient to rationally explain the different required energy during the reaction. 45Thus, the activation strain model (ASM) was employed to study the properties and behavior of reactants throughout the reaction process. 46The aim is to provide a comprehensive explanation for uncovering the mechanism of free radical scavenging reactions by antioxidant compounds.Besides studying the antioxidant capacity of syringol compound derivatives, drug-likeness predictions were conducted to evaluate their suitability as bioactive compounds upon consumption. 47[50]

Theoretical methods
Gaussian16 soware package was used for all density functional theory (DFT) calculations. 51The Minnesota function theory level, M06-2X, was utilized for every calculation performed during this study.This particular approach was chosen due to its ability to offer more precise and reliable geometric and energetic data compared to alternative DFT functions, notably the B3LYP function. 52,53The Pople basis set with diffusion and polarization functions, 6-311++G(d,p), was employed in this study due to its common application among researchers studying free radical scavenging activity. 54The density-based implicit solvation model (SMD) was performed in the calculations to mimic the physiological media for the free radical scavenging process. 55Water and pentyl ethanoate were used as representative media for body uids and lipid bilayer membranes of cells, respectively. 56,57n open-shell systems, unrestricted calculations have been conducted.Local minima and transition states (TS) are distinguished by the presence of imaginary frequency ( i f = 0 and i f = 1, respectively).Moreover, TS structures are also validated with intrinsic reaction coordinate (IRC) calculation to assure that pre-and post-complex are connected.Lastly, thermodynamic correction at 298.15 K was applied for all relative energies.Also, the standard state conversion into solution (eqn (1)) was included for the bimolecular reaction. 58Moreover, Okuno's correction (eqn (2)) is utilized likewise to consider entropy loss caused by the solvent cage effect. 30,59Those equations are described as follows.DG FV sol y DG The molar fraction of neutral and ionic species for each compound was determined by modifying the simple Henderson-Hasselbach equation into the following equations (eqn (3) and ( 4)) described below.For any compound whose experimental pK a values are still unknown, the parameter tting (PF) method was conducted.This method is proven to give similar results with the experimental data for phenolic compounds at any theoretical level. 60 In the thermodynamic and kinetic study, the QM-ORSA protocol was performed. 61In both the media, two main radical scavenging mechanisms were studied, namely, HT (Hydrogen Transfer) and RAF (Radical Adduct Formation).Only oxygen atoms in phenolic rings, benzylic, and allylic C sp 3 atoms, which could serve as hydrogen donors, were considered in the HT mechanism.Meanwhile, all unsaturated carbons were considered in modeling the RAF mechanism.Besides, Single Electron Transfer (SET) and Sequential Proton Loss Electron Transfer (SPLET) mechanisms were also included in polar media because of the presence of anionic species.On the other hand, these species could not be formed spontaneously in the lipid media and have been proven by prior work. 37Thus, only neutral species of syringol derivatives were studied in the lipid environment.
Only sites undergoing exergonic reaction (DG < 0) in each compound were followed by kinetic calculation.[64][65][66][67] where T is the temperature in Kelvin, R is the gas constant in kcal K −1 mol −1 , h is Planck's constant, k B is Boltzmann's distribution constant, and DG ‡ is the activation energy.9][70][71] For any scavenging mechanism that involved electron transfer, such as SET and SPLET, Marcus theory (eqn (6)) was performed to nd their activation energy ðDG ‡ ET Þ. 72,73 The value of reorganization energy (l) could be approximated using DG ET and DE ET (eqn (7)) values, where DG ET is the Gibbs energy of reaction, whereas DE ET describes the energy difference between the reactant and vertical product.Diffusion control was considered for any rate constant that exceeds 10 8 M −1 s −1 . 545][76][77] All equations that have been mentioned before are described as follows.
where k D is the steady state Smoluchowski rate constant, R AB is the reaction distance, D AB is the total mutual diffusion coefficient of free radical (A) and syringol derivatives (B), N A is Avogadro's number, h is the viscosity of the solvent, and r is the radius of the solute.For HT and RAF mechanisms, the reaction distance was measured from two heavy atoms in transition states, which were involved in the transfer process.Meanwhile, the reaction distance in SET and SPLET mechanisms was simply the total radius of solutes.The IRCs obtained were then used in activation strain analysis.Each potential energy (DE(z)) at the reaction coordinate point z was decomposed into strain destabilization energy (DE strain (z)) and interaction energy (DE int (z)), as detailed in eqn (11).The interaction energy was obtained from the difference between the complex structure energy and strain energy relative to the reactant complex. 78 Strain energy is the energy required to deform the reactant structure to adopt the structure present at each reaction coordinate point.Therefore, this energy value is always positive.According to eqn (12), strain energy was obtained by summing the relative single-point energies of the syringol derivative structure (DE strain syd (z)) and the radical (DE strain rad (z)) at each point on the IRC diagram. 79When performing single-point calculations, the syringol derivative structure was considered to have singlet spin multiplicity, while the HOOc radical has doublet spin multiplicity.
The three types of obtained energy are represented in a 2D diagram, with the distance of the key bond in the reaction as the abscissa.In the HAT mechanism, these are the broken bonds of O/H or C/H, while the formation of the C/O bond was observed in the RAF mechanism.

Acid-base equilibria
Based on the prediction results, the parameter tting method demonstrated reasonably good accuracy, although it is slightly lower than the actual values. 80The predictions also indicate that syringol and its derivatives are classied as weak acids.This can be related to the formation of intramolecular hydrogen bond between 1-OH with O-methoxy, which increases the deprotonation energy. 81,82At the physiological pH, syringol derivatives are highly available in a neutral form.As shown in Table 1, all of them exist only in the range from 0.23 to 0.48% as anions.Nonetheless, this amount is still considered in the subsequent analysis because it tends to react within the diffusion rate limits. 83,84ermodynamic study From the data in Table 2, all sites in the HAT mechanism are exergonic.The presence of substituents at the para (C-4) position appears to enhance the spontaneity of the reaction.Similar results have been found in syringol derivatives reported in other studies, such as canolol, 29 sinapinic acid, 31 and syringic acid. 32ydrogen abstraction from the phenolic sites reveals that the stability of the formed radical compounds can be ranked as follows: HPns > HPs z HAs > Hs.The position of the double bond in the propanoid group plays a crucial role to enhance product stability as the size of electron delocalization contributes to lowering the free energy of the reaction.
For carbon-centered radicals, the stability order is HAs > HPs > HPns.Unlike phenolic radical stability, an increased number of groups involved in resonance effects is an important feature to make the reaction more spontaneous.This is particularly noticeable in the HAs, which has both allylic and benzylic groups, allowing for greater electron delocalization.In contrast, the carbon radical in HPns is delocalized by only one group.
In polar solvents, antioxidant compounds may scavenge free radical activity through electron transfer mechanisms, specically SPLET and SET.However, we found that SPLET mechanism showed thermodynamic viability, whereas the SET reaction did not.These results differ slightly from other lignin derivative groups, such as guaiacol, where both mechanisms are less favored. 20Among the four studied compounds, electron transfer reaction in HPs is the most spontaneous.It appeared that the electron-donating power at the C-4 position inuences this reaction.
In contrast to the RAF mechanism, it is difficult for all syringol derivatives to form products in all studied solvents.The loss of aromaticity of the adduct formed leads to unfavorable scavenging activity by the phenolic ring.This effect was not observed in the propanoid group, where several reaction sites showed various ranges of viability.To gain an in-depth explanation of the ring's impact on the spontaneity of RAF mechanism, activation strain analysis was conducted and discussed in later section.
The position of the double bond in the propanoid group demonstrated a signicant inuence on the spontaneity of the Paper RSC Advances RAF mechanism.For instance, in HAs, the isolated sp 2 carbon atom separated from the resonance system leads to a decrease in the stability of the radical products.This results in a negligible impact of solvent polarity and charge on enhancing the delocalization effects, making an insignicant trend of spontaneity.

Kinetic study
Based on thermodynamic studies, it can be observed that the hydroperoxyl radical scavenging reaction of the syringol derivative occurred more spontaneously through the HT mechanism in all the solvents.However, kinetic analysis must still be carried out to quantify the contribution of each mechanism based on kinetic parameters, namely, activation energy and reaction rate constant.The transition structures at the HT and RAF sites, which are exergonic, are sought, and their activation energies are determined.The processed kinetic data has been summarized and can be seen in Table 3.
Consistent with the thermodynamic data, substituents also enhance the reactivity of syringol derivatives in all available reaction pathways.Apart from being examined through thermochemical analysis, this trend can also be characterized by the topological structure of transition states based on Hammond's postulate.For HT in the 1-OH mechanism, a higher order of the phenolic OH bond depicts an "earlier" transition state (TS), correlating with the lower activation energy.
However, compared to the propanoid group, this relationship could not be applied.As seen in Fig. 2, all carbon-centered reaction sites in all the derivatives have higher activation energy despite possessing an "earlier" TS structure.This is related to the formation of the hydrogen bond between the O-methoxy and the HOOc radical in the 1-OH transition states, which assists in stabilizing the structure, thus reducing the energy barrier. 85Therefore, the hydrogen transfer by 1-OH proceeds faster, despite its "later" TS structure.
The "earlier" TS characteristics also show a fairly good correlation with increased reactivity in the RAF mechanism.HPns, which has a smaller C-O bond order compared to HAs, possess an 'earlier' TS structure, resulting in a higher reaction rate.
However, despite having a higher activation energy, HPns actually exhibits a better radical scavenging rate through the RAF site compared to canolol in all solvents. 29It seems that the difference in the quantum tunneling effect on these two compounds underlies this phenomenon.
The polar environment also inuences the energy required for each site to react via the HT and RAF mechanisms.In an aqueous environment, all energy barriers for each reaction site are slightly lower than in a pentyl ethanoate environment.Decreasing energy barrier can also be observed from the 'earlier' TS at all sites in an aqueous solvent.
Among the two available forms, the anion species is the most reactive form in scavenging hydroperoxyl radicals in all the existing mechanism pathways.This species has lower energy barriers by 11 kcal mol −1 at all the sites and reacts within the diffusion limit range.However, it may not necessarily be dominant due to its availability, which affects its reactivity.
Table 3 Gibbs free energy of activation (DG ‡ ) in kcal mol −1 at 298.15 K, quantum tunneling coefficient (k), and realistic rate constant (k app ) for viable reaction pathways.HA: neutral species, A − : anionic species, PE: pentyl ethanoate solution, W: aqueous solution a Mechanism, site To obtain the overall rate constant (k overall ) for each reaction site in all the pathways in a non-polar solvent, the values of k app are summed following eqn (13).Meanwhile, in an aqueous solvent, these values are rst multiplied by the molar fraction of each syringol derivative species at physiological pH to obtain the molar fraction-weighted rate constant (k f ), which is then summed to obtain the overall rate constant in hydroperoxyl radical scavenging, as indicated in eqn ( 14)- (16).
In the case of Hs, its scavenging ability in a non-polar environment is solely inuenced by the HT mechanism at the 1-OH site because there are no other reaction pathways available.For the other three compounds that possess multiple available reaction pathways, the branching ratio (G i ) calculations were carried out to determine the inuence of each pathway on the scavenging reactivity using the equation below.
These values indicate that hydrogen transfer from the 1-OH site remains the sole pathway affecting the total rate constant, while other site contributions are negligible.Thus, the scavenging process of the HOOc radical through hydrogen abstraction from the phenolic site is a key factor in designing derivatives of phenolic-based antioxidants in a non-polar environment.This argument is supported by similar research on amino pyridinol derivatives, 92 sinapinic acid, 31 ferulic acid, 91 eugenol, guaiacol, and vanillin, 37 which show similar ndings.
In an aqueous environment, the HT mechanism at the 1-OH site also remains as the dominant pathway for nearly all the studied syringol derivatives in scavenging hydroperoxyl radicals.This domination was applied only within the acidic to physiological pH range, as seen in Fig. 4.Under basic conditions, SPLET becomes the key pathway in providing hydroperoxyl radical scavenging reactivity for syringol derivatives due to the increased availability of anionic species within that pH range.
The reactivity of syringol derivatives was also observed within a specic pH range, which is 2-12.In an acidic environment (pH = 2-6), all derivatives exhibited good peroxyl radical scavenging capacity, as seen in Fig. 5. Within this range, the reactivity order of the studied syringol derivatives was the same as that at the physiological pH, with HPs being the most reactive one.However, in an alkaline environment, there was a signicant reactivity increase and change in the order among the four syringol derivatives.HPns becomes the most reactive derivative one Table 4 Realistic rate constant (k app ), molar fraction-weighted rate constant (k f ), pathways, Branching ratio (G (%)), and overall rate constant (k overall ) of the studied syringol derivatives.HA: neutral species, A − : anionic species, PE: pentyl ethanoate solution, W: aqueous solution Mechanism, site among the others at pH = 8-12.Meanwhile, HAs and HPs react with hydroperoxyl radicals at nearly the same rate.Only syringol consistently exhibited less reactive scavenging ability compared to the other three compounds in every observed pH range.

Elucidating HT and PCET pathways
In antioxidant mechanisms, the hydrogen atom can be abstracted simultaneously in the form of protons together with electron, a process known as PCET (proton-coupled electron transfer).Therefore, the process yields the same products as HT, making it difficult to distinguish between them.Several methods have been proposed to differentiate these two mechanisms, although such a determination is not a straightforward matter.
The charge (Q) and spin density of transferred hydrogen atom (ASD) analysis at the transition state (TS) were performed according to the method reported by Galano and are summarized in Table 5.In the molecular form, the spin population on the migrating hydrogen atom at all the reaction sites indicated small negative values, consistent with the HT process. 93On the other hand, hydrogen transfer by the CH site in anionic species had a positive population value much greater than that reported previously.In connection with the charge analysis, these species can proceed through the HT mechanism due to their smaller charge value.However, this analysis method did not provide a satisfactory explanation as stated by Galano.
For this reason, SOMO (Single Occupied Molecular Orbital) topology analysis was then performed as it provides a reasonable description of the PCET mechanism.The observation indicated that the orbital densities at all the sites form nodes aligned to the hydrogen abstraction pathway, implying the HT mechanism.However, the cisoid form of TS at 1-OH obtained from the experiment showed that this analysis is less signicant for deducing it as the HT mechanism because the previous study reported that the HT mechanism occurs at those sites possessing planar or transoid TS structure, which affects our ndings. 94,95herefore, the observation was continued using HOMO (Highest Occupied Molecular Orbital) shape as recommended by DiLabio. 96As seen in Fig. S5, † all sites exhibited overlapping interactions between the p-orbital of the phenolic ring and the lone pair electron of the radical.This enables the transfer of electrons through different orbitals, leading to the PCET mechanism.The presence of this interaction also resulted in the orbital density of the SOMO not being the same as the topology of the conventional PCET mechanism. 97owever, the presence of similar overlap orbitals occurring at the benzylic site made the determination of the PCET mechanism a little doubtful as previous research found that those sites undergo the HT mechanism. 93,98Also, if HOMO analysis is ignored, SOMO analysis on 1-OH site gives different results from the previous study. 69Hence, we still could not give a nal hydrogen transfer mechanism for the radical scavenging activity of syringol derivatives.

Decomposition energy analyses by activation strain model (ASM)
Based on the previous and current studies, the reactivity of the HT mechanism on benzylic C-H is consistently lower compared to the phenolic OH, even though the reaction is generally  thermodynamically favoured due to the high activation energy.However, there is no in-depth study to explain this phenomenon.Therefore, we attempt to use the activation strain model to quantify the factors inuencing this value. 78s observed in Fig. 6, the destabilization strain energy at 1-OH is signicantly higher during the reaction compared to 1 0 -CH in all solvents.This is because of the intramolecular hydrogen bond that suppresses the movement of the phenolic H-atom to form the transition structure, as depicted in Fig. 7.This does not occur at the 1 0 -CH site, allowing the bond to transform freely with lower energy into the TS geometry.
However, during the reaction, the formation of a hydrogen bond between the HOOc radical and the syringol derivatives at the 1-OH site becomes a key factor in increasing the interaction energy, thus lowering the activation energy.The absence of this interaction at the 1 0 -CH site resulted in insufficient interaction energy to stabilize the strain energy, leading to a high activation energy.
Analysis was also conducted on the RAF mechanism to determine the factors causing the aromatic sites to react nonspontaneously with the HOOc radical.For this purpose, the potential energy during the reaction at the C-1 and C-2 0 sites were decomposed and compared.As presented in Fig. 8, the deformation of the alkene structure requires lower energy compared to the aromatic ring.Besides the loss of p electron delocalization, the high deformation energy at the C-1 site is also due to the additional energy required to disrupt the relatively rigid planarity of the aromatic ring.
Surprisingly, the formation of radical adducts by aromatic sites has higher interaction energy.It appears that there are non-covalent interactions between the HOOc radical and the aromatic ring during the transition structure formation.The solvent polarity also seems to assist in enhancing these interactions, making the formation of transition structures in aqueous solvent more favorable.
However, the loss of ring aromaticity aer the TS formation resulted in a decrease in the interaction energy, which is insufficient to stabilize the increasingly higher strain energy.It is also shown that the interaction energy plot also attens,  leading to an intersection with the interaction energy plot of C-2 0 .This resulted in the formation of an unstable product that can be readily dissociated back into reactants.

Biocompatibility prediction
Although all the studied compounds have been proven to be effective antioxidants, the physicochemical and toxicity parameters need to be considered to determine whether these four compounds are suitable for consumption.Therefore, the SwissADME program was employed to calculate the physicochemical and pharmacokinetic properties that are assessed in the ADMET analysis. 99The parameters used refer to the Lipinski, 100 Ghose, 101 and Veber 102 rules, as detailed in Table 6.Meanwhile, the LD 50 values and mutagenic potential (M) are used to assess the toxicity of syringol derivatives, calculated using the T.E.S.T (Toxicity Estimation Soware Tool) program. 103,104ll the obtained values are then summarized in Table 7. From the ADME analysis, all syringol derivatives have met the necessary criteria.These syringol compounds have sufficient molecular masses to easily pass through the cell membranes.Although the molecular mass of Hs slightly violates Ghose's rule (MW = 154.16), it can be tolerated, as many commercial drugs, such as isoniazid, allopurinol, and metformin, have molecular masses < 150 g mol −1 . 101he cell permeability of syringol derivatives can also be predicted from the values of lipophilicity, polar surface area, and molar refractivity.All syringol derivatives have log P values in the range of −0.4 to 5.0, indicating good solubility in both   polar and non-polar environments.The number of hydrogen donors and acceptors, as well as the total number of elements in the compounds, also affect these properties.The interesting part is the moderately lipophilic nature and the polar surface area of syringol derivatives (0.4-0.6; <79 Å 2 ), which allow them to pass the blood-brain barrier (BBB), as indicated by the BOILED-EGG model. 105This ability was also validated by making predictions using another program called LightBBB, which yielded similar results. 106The ability to pass the membrane is known as a requirement for the development of neuroprotective compounds. 107o gain an understanding of the potential hazards of syringol derivatives on the body, the LD 50 and M values of commercial antioxidants (ascorbic acid and melatonin) were used as comparative data.Substituted syringols are known to have higher LD 50 values than melatonin, while Hs has an LD 50 value of 549.58 mg kg −1 , which is signicantly below the reference values.Other studies have also reported that syringol shows toxicity effects in vitro and in vivo in the concentration range of 0.5-2 ppm. 108Additionally, Hs exhibits potential to induce mutations in the body (M = 0.53), which is not observed in other syringol derivatives.Considering all these factors, apart from their potential as safe antioxidants, compounds like HAs, HPns, and HPs can be further investigated as specic neuroprotective agents.

Conclusions
Kinetic and thermodynamic studies using the M06-2X/6-311++G(d,p) level of theory of compounds Hs, HAs, HPns, and HPs have demonstrated good antioxidant capacity in both polar and non-polar environments.Among these four compounds, HPs is the most effective as HOOc radical scavenger with k overall values of 1.75 × 10 8 M −1 s −1 in water and 1.13 × 10 8 M −1 s −1 in pentyl ethanoate.In pentyl ethanoate, hydrogen transfer at the 1-OH site is the most contributing pathway in all syringol derivatives.Observations of the transition structure at the 1-OH site showed the presence of intermolecular hydrogen bonding, which inuences the low activation energy.This interaction is quantied using the Activation Strain Model (ASM), which has been successfully employed for the rst time in antioxidant activity research.Furthermore, the ASM analysis can also explain the disfavor of the RAF mechanism by the aromatic ring, as evidenced by the decrease in the interaction energy aer the formation of the transition structure.Hydrogen transfer mechanism in syringol derivatives is still a challenge to ensure whether through HT or PCET reactions.Under physiological conditions, HT competes with SPLET reactions (43.4% vs. 56.6%) in Hs, with SPLET being the dominant pathway.Ionization effects and pK a values contribute to the high radical scavenging capacity of HPns in basic environments.Based on physiochemistry and toxicity parameters prediction, the four compounds exhibited good bioavailability and low toxicity.Only HAs, HPns, and HPs had the potential to be used as bioactive compounds in regard of health.The BBB-permeable properties of these three compounds also provide insights for further research in the therapy of neurodegenerative diseases.Although Hs is toxic to humans, its antioxidant abilities may nd applications in other industries.

Fig. 2
Fig. 2 Transition state structures of the HT mechanism at the O-phenolic site in each syringol derivatives, along with bond lengths (in black), bond angles (in black), and Wiberg bond orders (in blue) in water and pentyl ethanoate (in bracket).

Fig. 3 Fig. 4
Fig. 3 Reactivity power of the studied syringol and other known antioxidants in (a) non-polar media and (b) polar media relative to Trolox.

Fig. 5
Fig. 5 pH influences in the overall rate constant of syringol derivatives.

Fig. 6
Fig. 6 Activation strain model of the HT mechanism of 1-OH (black) and 1 0 -CH (red) in (a) water and (b) pentyl ethanoate using HAs as the model.

Fig. 7
Fig. 7 Hydrogen bond in (a) the reactant complex and (b) the transition structure in HT mechanism of HAs.Bond length in the studied media is shown as water (pentyl ethanoate).

Fig. 8
Fig. 8 Activation strain model of the RAF mechanism between C-2 0 (black) and C-1 (red) in (a) water and (b) pentyl ethanoate using HPns as the model.

Table 1
pK a values and molar fraction of syringol derivatives at physiological pH.HA: neutral species, A − : anionic species a M f A − (%) a Obtained from ref. 80.Table 2 Gibbs free energy of the reaction (DG) in kcal mol −1 at 298.15 K for the modeled reaction pathways.HA: neutral species, A − : anionic species, PE: pentyl ethanoate solution, W: aqueous solution © 2024 The Author(s).Published by the Royal Society of Chemistry RSC Adv., 2024, 14, 6310-6323 | 6313

Table 5
Atomic spin densities (ASD) and charge (Q) of migrating H atom from natural orbital population analyses.HA: neutral species, A − : anionic species, PE: pentyl ethanoate solution, W: aqueous solution

Table 6
Physicochemical analysis parameters and their values to satisfy Lipinski's, Ghose's, and Veber's criteria

Table 7
Physicochemical parameter of investigated syringol derivatives and relevant drugs